Elastic medium confined in a column versus the Janssen experiment. 
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We compute the stresses in an elastic medium confined in a vertical column, when the material 
is at the Coulomb threshold everywhere at the walls. Simulations are performed in 2 dimensions 
using a spring lattice, and in 3 dimensions, using Finite Element Method. The results are compared 
to the Janssen model and to experimental results for a granular material. The necessity to consider 
elastic anisotropy to render qualitatively the experimental findings is discussed. 

PACS numbers: 

I. INTRODUCTION 

The mechanical status of granular matter is presently one of the most open and debated issues pj. This state 
of matter exhibits many unusual mechanical and rheological properties such as stress induced organization at the 
microscopic or at the mesoscopic Q level which may yield macroscopic effects such as stress induced anisotropy 
0,0]. This issue sets fundamental questions relevant to the understanding of many other systems exhibiting jamming 
such as dense colloids or more generally soft glassy materials 0,0- For practical applications, the quasistatic rheology 
of granular assemblies is described using a phenomenological approach, based on an elasto-plastic modelling of stress- 
strain relations 8J. So far, there is no consensus on how to express correctly the macroscopic constitutive relations 
solely out of microscopic considerations and under various boundary conditions or loading histories. This very basic 
issue was illustrated in a recent debate on how to understand the stress distribution below a sand pile and especially 
how to account for the dependence on preparation protocols |9| . A new mechanical approach was proposed based on 
the concept of " fragile matter" and force chains propagation modelling . But recent experiments have dismissed 
this approach and evidenced results more consistent with the traditional framework of general elasticity [Til ] . 

In this paper we focus on the predictions for stresses measurements at the bottom for an elastic material confined 
in a rigid cylinder. When the column is filled with granular material it corresponds to the classical Janssen's problem 
[l^. Recently, this issue has received a lot of attention either experimentally [131 IbH ITEl ITil fill fl8L Il9l ] or numerically 
|2d l2l| . the confined material being either pushed or pulled down. Surprisingly, so far to our knowledge, there are 
very few systematic comparison or even direct relation with the outcome of standard elasticity in the same situation 
of confinement. Note that experimentally, it was found that essentially elastic materials like gels, may display a 
Janssen stress saturation effects that could well predict the onset of self-collapsing under gravity |22|. In a recent 
paper [23}, Evesque and de Gennes proposed a model for the slow filling of an elastic medium modelling a granular 
packing. As the material is poured in the column, displacements of the material at the bottom are induced by the 
weight of the new material added so that it mobilizes friction forces. Within the assumption of a minimal anchorage 
length, it leads to partial and inhomogeneous mobilization of friction at the walls: friction is fully mobilized at the 
bottom, and partially in the upper part of the column. A Janssen like pressure profile can be derived in the case 
where the saturation length A is high compared to the column radius R; this last condition is actually restrictive and 
inappropriate for usual experimental cases |l3lll6|. 

We propose to study in detail the elastic predictions and to compare them to the Janssen model and to experimental 
results for a granular material in the same situation i.e. when the material is at the Coulomb threshold everywhere 
at the walls. The main comparison features with the experimental results have already been presented in |l6|: here, 
we detail much more the elastic predictions. 

Two kinds of situations are considered. First, the mass at the bottom of the column is measured as a function of 
the material filling mass. Second, similar measurements are produced with an overweight on the top of the material. 
We recently performed the corresponding experiments for a granular material and it was shown that one could 
obtain quite reproducible data provided a good control of the packing fraction homogeneity and a polarization of all 
the friction forces at the walls in the upwards direction. Our measurements confirmed, for the first time, the general 
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FIG. 1: Apparent mass M a vs. filling mass Mfm for a loose packing (y = 59%) of slightly polydisperse 1.5 mm glass beads in 
a 38 mm diameter steel column of friction coefficient p s — 0.25 without an overweight (squares), and with a 80.5 g overweight 
(triangles, line) on top of the granular material. The results of both experiments are compared to the Janssen model prediction 
(dotted lines). 

validity of Janssen's saturation curve. They also evidenced an overshoot effect of spectacular amplitude induced by a 
top mass equal to the saturation mass. These experimental results are actually strong tests for any theory of granular 
matter. 



II. EXPERIMENTAL RESULTS AND JANSSEN MODEL 

We first summarize our recent experimental results ^|| and compare them to the Janssen model predictions. In 
|l6j . a slightly polydisperse assembly of 1.5 mm glass beads was poured at controlled packing fraction in steel columns 
of various friction coefficient, and of diameter varying between 38 mm and 80 mm. The aim of the experimental 
procedure is to achieve the Coulomb threshold everywhere at the walls when the apparent mass M a is measured; see 
|l6j for details on the procedure. The apparent mass M a measured at the bottom of the column was plotted as a 
function of the filling mass Mfm of the material. The typical results obtained when Mfm is varied are shown on 
Fig. ^ The apparent mass M a saturates exponentially with Mfm. When an overweight equal to the saturation mass 
M sat is added on top of the granular material, M a increases with Mfm , up to a maximum M max , which is about 20% 
higher than M sat , then decreases slowly towards the saturation mass M sat . 

The simple model which captures the physics of this saturation phenomenon was provided in 1895 by Janssen [T^ | . 
This model is based on the equilibrium of a granular slice taken at the onset of sliding everywhere at the walls; we 
attempted to realize as best as possible this last condition in our experiment [T^ |. In cylindrical coordinates with 
origin at the top surface and the cylinder axis being the z axis, the relation, at the slipping onset, between the shear 
stress a rz and the horizontal stress a rr at the walls is 



a rz (r=R, z) = fj, s a rr (r=R, z) 



(1) 



where fi s is the Coulomb static friction coefficient between the grains and the walls. It results in a relation between 
the filling mass Mfm and the apparent mass at the bottom M a of the form: 



M =M sat (l-exp(-^)) 



M,, 



with 



M,at = 



pnR 3 
2KjT s 



(2) 



(3) 



where p is the mass density of the granular material, and K is the Janssen parameter rendering the average horizontal 
redirection of vertical stresses: 



Ka z 



(4) 
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From a mechanical point of view, a major simplification of this model comes from the assumption that the redirection 
parameter K would stay constant along the vertical direction. But on the other hand, it provides a clear and simple 
physical explanation for the existence of an effective screening length A = R/2K/i s above which the mass weighted 
at the bottom saturates. 

In [l6j . several saturation profiles were measured for various packing fractions, columns sizes and friction coefficients 
between the grains and the walls. When the apparent mass rescaled by the saturation mass is plotted as a function 
of the filling mass also rescaled by the saturation mass, we obtain a universal reseating of all data on a curve which is 
precisely the one predicted by Janssen (Fig. QJ: M a /M sat = f(Mfut/M sa t), with f(x) = 1 — exp(— x). The rescaling 
with radius R and friction coefficient /i s was also checked, and good agreement with the Janssen model rescaling was 
found. The Janssen constant K was found to depend on packing fraction v and an effective relation was derived: 
AK/K ~ 5Av/v. On the other hand, when a top mass equal to the saturation mass is added on the top of the 
granular material, the apparent mass M a displays a maximum M max 20% higher than M sat , whereas the Janssen 
model predicts M a = M sat whatever the filling mass Mfui is. Therefore, this overshoot goes beyond the possibilities 
of Janssen's model which seems adapted to a unique configuration. 

In the next section, we study in detail the predictions of isotropic homogeneous elasticity. 



III. ELASTICITY: THEORY AND SIMULATION METHODS 



A. Theory 

We first recall the general framework of homogeneous isotropic elasticity, and then predict the behavior of an elastic 
material confined in a column. 

The elastic theory gives, in the limit of small deformations, a linear relation between the stress tensor components 
aij and the strain tensor components e^ . For an isotropic elastic material, we get 

Eeij = (1 + v v )a l j — VpSijdkk (5) 

where E is the Young modulus, and v v the Poisson ratio which takes its value between —1 and 1/2 in 3D, and between 
-1 and 1 in 2D. 

In a uniaxial homogeneous compression experiment (Fig. |3J), where a zz — —p is imposed everywhere, the other 
stress tensor components being null, we get 

e« - -P/E (6) 

everywhere and 

£xx — £yy — Vp€ zz (7) 

which signifies that the material expands in the transverse direction. 

The Young modulus E is thus characteristic of the material's stiffness; a cylinder of length I and section S has stiffness 
k = ES/l in the axial direction. The Poisson ratio is linked to the material compressibility: the volume variation is: 
SV/V = —(1 — 2v p )p/E. Therefore, an incompressible material has Poisson ratio 1/2 in 3D (1 in 2D). 

If we now confine an clastic medium of Young modulus E and Poisson ratio v p , in a rigid cylinder of radius R, no 
more radial displacement at the walls is allowed: u r {r = R) — 0. Stresses and displacements can actually be calculated 
in the limit of high depths z under the assumption that they then should be independent of z. The boundary conditions 
we impose are the Coulomb condition everywhere at the walls 

cr rz (r=R) = fj, s a rr (r=R) (8) 

and infinitely rigid walls i.e. u r (r=R) = 0. The stress tensor components are then 

a rz (r,z) = —pgr (9) 

(10) 



l-Vp 

(1 - v p )pgR 



gff(r,s) = ~* - (11) 
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FIG. 2: Sketch of uniaxial compression of a free elastic material (top) and of a confined elastic material (bottom). 



And the asymptotic displacements are 

u z (r, z) = ^—S-pqr 2 -pqRz + Un (12) 

u r (r, z) = ug(r, z) = (13) 

This can be checked by injecting this solution in the stress-strain relation (JjJJ and internal equilibrium relation 

di<7 i: j = -pg (14) 

Thus, we obtain a Janssen's like redirection phenomenon due to a Poisson's ratio effect with a local Janssen's 
parameter K e i = a rr (r, z)/a zz (r, z), being for large depths: 

K el = (15) 

t Vp 

At 2D, we obtain the same saturation and stress redirection phenomena with K e i = v p . 

For a free elastic medium, the Poisson ratio effect is a transverse dilatation; for a confined elastic material, the 
Poisson ratio effect is a transverse redirection of stresses (Fig. (2J. 

In the following, K e i will design the elastic stress redirection constant, whereas K is devoted to design the stress 
redirection constant in the Janssen framework. 

In this section, we obtained the asymptotic values of stresses and displacements. For the vertical stress, the limit 
is similar to the Janssen asymptotic vertical stress if one identifies K an K e \. But we cannot say anything from this 
calculation about the whole pressure profile. We thus need to perform numerical computation. In the next section, 
we present the numerical methods we employed to simulate Janssen experiments for an elastic column. 



B. Numerical methods 



Two different methods were employed: we first computed the stresses in 2D with a spring lattice. We also computed 
the stresses in 3D, using Finite Element Method thanks to CAST3M [24|. 



1. 2D: spring lattice 



The 2D computations are performed in order to provide elastic predictions for a direct comparison with 2D simu- 
lations of granular materials. 



FIG. 3: Sketch of the discrete 2D elastic medium. 



In order to avoid any confusion, let us recall that we are not trying here to give a microscopic model for a granular 
material, but we use a discrete system behaving like an effective elastic medium in the continuous limit in order to 
perform simple numerical simulations. 

However, this system can describe the most simple of granular materials: a 2D hexagonal piling of frictionless disks 
(with non-hertzian contacts so that contact elasticity is linear). Note also that systems of springs have been recently 
studied by Goldenberg and Goldhirsch j2j|; they showed that large forces inhomogeneities (i.e. like forces chains) can 
be found at the discrete scale in these systems, which are elastic in the continuous limit. 

In order to simulate a 2D elastic medium, we put point masses m on an hexagonal lattice of link size a (Fig. 
|3J). Every particle is linked to her 6 neighbors with identical springs of stiffness k and length a at rest. Therefore, 
the potential interaction energy between particles placed at and Xi (such that ) at rest, submitted to 

infinitesimal displacements Ui and is: 

E P = \ k {^/ ( Xi + Vi ~ U i)( X i + v i ~ u i) ~ °) ( 16 ) 

By varying the stiffness k, we can only vary the young modulus E of the effective elastic medium. As we also need 
to vary the Poisson ratio v, an elastic torsion potential between neighbor springs separated by angle 6 (Fig. 0} is 
added: 

E p = h b cos 2 (17) 

This potential tries to maintain an angle 7r/3 between 2 neighbors springs if kb < 0. Thus the limit fcf, — > infty 
corresponds to a contractive medium of Poisson ratio -1, whereas the limit fc& — > oo corresponds to a incompressible 
medium. 

An elementary area A = a 2 v / 3/2 can be associated to each particle. The surface energy u> then reads: 

2 (\ 

w = ( g ^2 E p( s P rin 9 s ) + ^2 E p( an 9 l es) 

If we consider the continuum limit of this system, the displacement of a particle located at X4 is: 

Ui = {djUi)xj (19) 
The surface potential energy can be easily computed at second order in = (l/2)(diUj + djUi), and one obtains: 

/g / 3 1 \ 

u = T \4 {k + 2fc ' )(e ^ + (2yy) + 2 (fc " 6kb ^ e yy + {k + 6k ") e ly) ( 20 ) 

For an elastic medium, stress-strain linearity reads 




°y — Cijki^ki 



(21) 
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FIG. 4: Displacements of the points of an elementary cell of the 2D elastic medium. 



There are a priori 9 independent parameters in Cijki (but only 2 if it is an isotropic elastic medium). The surface 
energy then reads 

If we identify this energy to the one we obtained for the spring lattice, we see that our system, in the continuous 
limit, is an isotropic elastic medium of Young modulus E and Poisson ratio v v : 

1 fc — 6fc b 

^ = 3kTW b (23) 

3 k + 2k b v ' 

In order to solve the equilibrium problem of the lattice, we need to express the internal forces. They can be deduced 
from the surface energy l|22|) : to the usual elastic forces due to compression (or decompression) of springs, add forces 
which tends to drive the triangles angles to their equilibrium value (if k b < 0) or out of the 7r/3 value (if k b > 0). The 
force on A, due to the out of equilibrium angles of triangle ABC (Fig. 0} reads: 

3 

F A, = -7> k b{u Ai +R- w /3 i UB j + R 7r / 3 U Cj ) (25) 

Z 3 3 

where Rg ij is the rotation matrix of angle 9. 

For the numerical computation, we impose the balance of forces (gravity, elastic forces, torsion forces) everywhere. 
At the bottom, we impose a null vertical displacement (rigid bottom), and either a perfectly stick (i.e. nullity 
of horizontal displacement) or perfectly slip (i.e. nullity of horizontal projection of forces) bottom. At the walls, 
we impose a null horizontal displacement (rigid wall), and we impose the Coulomb condition: the forces projected 
vertically are proportional to the forces projected horizontally with proportionality factor fi s . On top of the material, 
the forces projected horizontally are null (perfectly slip overweight), and the overweight is simulated by imposing the 
same vertical displacement for each particle on top (stiff overweight), which is close to the experimental situation. We 
thus obtain a linear system on the point displacements. We can vary stiffness k and k b in order to simulate elastic 
mediums of different Young modulus and Poisson ratio. Note that for varying v v between and 1, we need to vary 
k b between — k/6 and fc/6. We also vary the friction at the walls [i s . 



2. 3D: FEM 

In order to get the whole stress saturation curve, finite element numerical simulations [24[ were performed. The 
column is modelled as an isotropic elastic medium. We vary the friction at the walls /i s , the Young modulus E and 
the Poisson ratio v p . We imposed a rigid (nullity of vertical displacements), either perfectly stick (nullity of horizontal 
displacements) or perfectly slip (nullity of horizontal stresses) bottom. We found no appreciable difference between 
these two previous cases. The condition a rz = [i s a rr is imposed everywhere at the walls. The cylinder is modelled as 
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a steel elastic medium. We verified that in all the simulations performed, there is no traction in the elastic medium, 
so that this can actually be a model for a granular material. 

For the simulations performed without overweight, the top surface is set free (no stress); the overweight is modelled 
as a perfectly slip (no horizontal stress) brass elastic medium. 

In order to impose the Coulomb condition at the walls, we first set F z = at each point of the mesh at the walls; 
we then obtain in these conditions the value F r exerted by the elastic medium on the walls. We then iterate in order 
to obtain F z = fi s F r where fi s is the friction coefficient at the walls: the vertical force imposed at step (i+1) is: 

F z (i + 1) = (1 - e) x F,(i) + e x fi s F r (i) (26) 

We choose e = 0.2. This procedure ensures convergence towards the Coulomb condition: if at step (i) the Coulomb 
condition F z (i) = /j, s F r (i) is fulfilled, then at step (i+1): F z (i+1) = (l~e)F z (i)+e/i s F r (i) = (l — e)n s F r (i)+efj, s F r (i) — 
fj, s F r (i) — F z (i). This boundary condition is the same as the one at step (i): this yields F r (i + 1) = F r (i), and thus 
F z (i + 1) = n s F r (i + 1). 

3. Remarks 

Note that in these simulations we imposed the Coulomb condition everywhere at the walls. This allows comparison 
with the Janssen model in the same situation. But, regarding the experimental results |Tfij . i) nothing really insures 
that our piling preparation is strictly isotropic and ii) in spite of the careful procedure, we are never absolutely sure 
that all the friction forces at the wall are actually mobilized upwards. Moreover, the modelling of the contacts may 
seem rudimentary, as elasticity of contact should be included. 

In these simulations, imposing dynamical friction at the walls or the static Coulomb threshold is actually the same: 
the material obeys the same equilibrium equations (if we consider a steady sliding at the walls), and the same condition 
at the walls, with just a change in the name (and the experimental value) of the friction coefficient. 

IV. SIMULATION RESULTS 

In this section, we present the results obtained from numerical computations at 2D and 3D. 

In the following, we vary mainly the friction coefficient, and the Poisson ratio. The 2D simulations are in arbitrary 
units. The 3D simulations were all performed, if no contrary mention, for an isotropic elastic medium of mass density 
p = 1.6 gem -3 (which corresponds to an assembly of glass beads of packing fraction D — 64%), of Young modulus 
100 MPa, in a steel cylinder of radius R — 4 cm, Poisson ratio 0.3, Young modulus 210 GPa, and thickness 3 mm. 
These data, which correspond to the display used in |Tfjj . will not be specified anymore in the following. Simulations 
will also be performed in order to study the influence of the variation of the cylinder radius and the Young modulus 
of the elastic medium. 



A. Simulations without overweight 

We first perform simulations similar to the original Janssen experiment: we plot the weight at the bottom as a 
function of the weight of elastic material in the column (Fig. [5j|. 

We see on Fig. \5\ for friction coefficient fi s = 0.5 and Poisson ratio v p = 0.45, that the apparent mass M a saturates 
exponentially with the filling mass. The data are perfectly fitted by the Janssen model for these parameters (Fig. 
|SJ), but the Janssen coefficient K extracted from the fit is 9% higher than the elastic stress redirection constant 
K e i — f p /(l — Vp) — 0.82. This is a priori unexpected as the elastic saturation pressure and the Janssen one should 
be identical with K = K e i. 

1. Effect of the bottom 

In order to understand this feature, we study the whole mean vertical pressure profile in the same column (Fig. ■ 
Regarding mean vertical pressure, we see that the asymptotic value is now the expected theoretical asymptotic 
value, and the Janssen curve with K = K e i gives a good though not perfect fit of the profile. The pressure saturates 
at high depths but decreases suddenly near the bottom; this is actually the value on the bottom we measure in a 
Janssen experiment. This feature explains why the saturation mass is lower than the expected one on Fig. [3] The 
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Weight of the elastic medium in the colulmn (N) 



FIG. 5: 3D simulation of a Janssen experiment (squares) for an elastic material of Poisson ratio v v — 0.45; the friction at the 
walls is fis = 0.5. The data are fitted by a Janssen curve of coefficient K = 0.89 (line). The Janssen curve for K corresponding 
to the elastic stress redirection constant K e i = v p /(l — v p ) — 0.82 is also displayed (dotted line). 
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FIG. 6: Mean pressure profile in a simulated elastic material of height 31 cm and Poisson ratio v v = 0.45 (squares); the friction 
at the walls is (is = 0.5. Depth z = cm corresponds to the top of the column. Depth z — 31 cm corresponds to the bottom. We 
display the integral of vertical stresses F z at height z. The data are compared to a Janssen curve of coefficient K — K e i = 0.82 
corresponding to the elastic stress redirection constant (dotted line). 



reason for this change of pressure near the bottom is that the asymptotic vertical displacement is parabolic whereas 
the bottom is rigid and imposes a fiat displacement. 

It is thus important to note that the usual Janssen experiment, in which one measure is made for one given height, 
is not equivalent to measuring a pressure profile, and results in a lower saturation stress (i.e. higher Janssen's constant 
K) than the pressure profile. For more clarity on this last point, we illustrate this difference on Fig. \7\ 
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FIG. 7: Sketch of comparison between a pressure profile and measures at the bottom for an elastic medium. 
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FIG. 8: Study of the rescaling with the Poisson ratio v v for 3D simulations of an elastic material. We plot M a x fj, 3 v p /(l — v v ) 
vs. Mfui x fJ, s Vp/(l—i/p) for Up — 0.26 (squares), v v — 0.35 (open circles), v v — 0.4 (triangles), v v — 0.45 (open down triangles) 
and Vp — 0.49 (stars); the friction at the walls is fi 3 — 0.5. 
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FIG. 9: Janssen constant K extracted form 2D (squares) and 3D (triangles) simulations for various Poisson ratios. In the 2D 
simulations, the friction at the walls is fi s = 1.0; in the 3D simulations, the friction at the walls is /i s = 0.5. K is plotted vs. 
the elastic stress redirection constant K e i = v p in 2D, K e i = v v j(\ — v v ) in 3D. 



In the following, we study the rescaling law of simulated Janssen experiments with friction coefficient /i s at the walls 
and the Poisson ratio v v . The data are compared to the Janssen model predictions by plotting M a x fi s x v p /(l — v p ) — 
f(Mftii x fi s x v p /{l - v p )) for different ^ s and v v (since K e i = v p /{1 - v p ) and M sat (xl/(K e ifi s )). 



2. Effect of the Poisson ratio 



On Fig. |H1 we study the rescaling with the Poisson ratio v v for 3D simulations. The rescaling law is rather good, 
the differences may not be observable experimentally. On Fig. [5]we plot the Janssen coefficient K extracted from the 
fit of data in 2D and 3D versus the elastic stress redirection constants K e i {y v in 2D, v v j(\ — v v ) in 3D). 

We observe that for a given friction coefficient, K hardly depends on v p in 2D: K variation is 1% for K e i varying 
from 0.33 to 0.77. In 3D, K/K e [ increases roughly linearly with K e i for K e i varying from 0.25 to 0.96; however, K 
variation is less than 10% in this range. We remark that K > K e f. we explained it in the preceding section as a 
consequence of the presence of the rigid bottom. 
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FIG. 10: Study of the rescaling with the friction coefficient fi s at the walls for 2D simulations of an elastic material. We plot 
M a x M» vs - Mfai x /is for fi s = 0.4 (squares), fi s — 0.6 (open circles), fi s = 0.8 (triangles) and Us = 1.0 (open down triangles); 
the Poisson ratio is v v = 0.77. 
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FIG. 11: Study of the rescaling with the friction coefficient [i s at the walls for 3D simulations of an elastic material. We plot 
M a x fj, s v P /(l — Vp) vs. Mfui x fj, s Vp/(l — Vp) for /i s = 0.1 (open down triangles) fi s = 0.25 (triangles), fi 3 = 0.5 (open circles) 
and /i s = 0.8 (squares); the Poisson ratio is v v — 0.45. 



3. Effect of friction at the walls 

On Fig. EJand Fig. ^2 we study the rescaling with the friction coefficient fj, s at the walls respectively for 2D and 
3D simulations. We now see that the proposed rescaling law is not good: the saturation is more abrupt as the friction 
coefficient is higher. 

On Fig. El w e plot the Janssen coefficient K extracted from the saturation mass value in 2D and 3D versus fi s . 
We now observe that the Janssen constant (and thus the saturation mass) depends strongly on the friction coefficient 
at the walls: K increase is 12% in 2D for fi s varying between 0.1 and 1.0, and 20% in 3D for fi s varying from 0.1 to 
0.8. K depends roughly quadratically on fj, s . Moreover, the Janssen constant seems to tend towards the elastic stress 
redirection constant K e i when fj, s goes to 0. 

In order to understand these features, we plot on Fig. 1131 the mean vertical pressure F z versus the depth, and its 
rescaling with /i s . 

We see that, regarding the pressure profile, the Janssen rescaling law is correct for the asymptotic value, but the 
profiles are slightly different: for low friction at the walls fi s as 0.25, the profile is perfectly fitted by the Janssen law; 
for higher friction, the profile is sharper and saturates abruptly. The differences for simulated Janssen experiments 
observed on Fig. El are now identified as a consequence of the presence of a rigid bottom: we actually see on Fig. El 
that the decrease of pressure near the bottom is more important when the friction is higher i.e. the saturation mass 
is lower (and the effective Janssen constant is higher). 
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FIG. 12: Janssen constant if extracted form 2D (squares) and 3D (triangles) simulations for various friction coefficient pi s . In 
the 2D simulations, the Poisson ratio is v v — 0.77; in the 2D simulations, the Poisson ratio is v v — 0.45. K is rescaled by the 
elastic stress redirection constant K B i = 0.77 in 2D, K e i = 0.82 in 3D. 
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FIG. 13: Integral of vertical stresses F z vs. depth z, for elastic materials of Poisson ratio v v — 0.45, and friction at the 
walls fis = 0.25 (triangles), fi s = 0.5 (open circles), et fi s = 0.8 (squares). We plot F z /F sa t vs. 2K e ijj, 3 z/R where F sa t = 
pgirR 3 /(2K e ifi s ) is the theoretical saturation value for F z . z/R — corresponds to the top of the material. We also plot the 
Janssen model prediction with K = K e i — v v /{\ — v p ) (line). 

We now understand all these features: the parabolic part of asymptotic displacements is negligible for low friction; 
in this case, the flat displacement imposed by the rigid bottom matches the asymptotic displacement, and the pressure 
at the bottom is the saturation pressure: we thus obtain K = K e i . As the friction is increased (and the material leans 
on the walls), the parabolic part of asymptotic displacement becomes more important and the influence of the rigid 
bottom is to decrease pressure; we thus observe an increase in the effective Janssen constant K with fj, s . 



4- A unique parameter 



In the Janssen picture, the Janssen coefficient K and the friction at the walls /i s are 
the relevant parameter is K/i s . As regards the saturation mass M sa t, it remains true in 
I14l that Kfi s , extracted from the saturation mass value, is a function of K e i\x s alone (i.e. 
replotted on a single universal curve); we find a quadratic dependence of K\x s on K e ifi s 
This explains the dependence found before on v v at fixed fi s or on ii s at fixed v v . 

However, it is not true anymore for the whole Janssen profile: we see on Fig. II 51 that there is not a universal curve 
M a /M sat = f(Mfui/M sat ): if most data may be fit by the Janssen curve (with f(x) = 1 — exp(— x), it is not true 
anymore for high (i.e. low M sat values); the data for high values of K[i s (=0.65 here) fall above the Janssen curve. 



not independent parameters: 
elasticity, and we see on Fig. 
data on Fig. I§l and 112*1 can be 

K^ s « K el n s + 0.29(K P jn s ) 2 . 
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FIG. 14: Janssen constant K extracted form 3D (triangles) simulations for various friction coefficient fi a at constant v v = 0.45 
(empty triangles) and various Poisson ratios v v at constant fi 3 = 0.5 (squares); see Fig. I%1 and fTTI for v v and \i a values. The 
dotted line is the y — x line; the line is a polynomial fit with Kfi B — K e ifi a + 0.29(K e i/i s ) 2 . 
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FIG. 15: Ma/Msat vs. Mfai/Msat for 3D simulations for various values of v v and fj, s ; the triangles are for v v = 0.45 and 
Us = 0.8. 



However, for most experimental conditions, K/i s is not as high, and this effect will not be observable. 



5. Effect of walls elasticity 

On Fig. 1161 the Young modulus is varied. 

We observe that the results are independent of the Young modulus E value, as long as it is less than 500 MPa; we 
see on Fig. ^|that the curves obtained for the simulation of a Janssen experiment for E = 1 MPa and E = 200 MPa 
can be perfectly superposed. For higher E, the weighted mass is higher and does not seem to saturate anymore; as 
an example, for E = 4 GPa on Fig. 1161 M a seems to increase indefinitely. 

Note that this Young modulus effect is present only because we take into account the walls elasticity; for rigid walls 
and bottom, there would not be any dependence on E. The values of E presented here have a meaning only for a 
particular cylinder (same E and same thickness). 
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FIG. 16: Simulation of a Janssen experiment for elastic materials of Young modulus E = 1 MPa (squares), E — 200 MPa (open 
triangles) et E — 4 GPa (circles), of Poisson ratio v v — 0.45, in a steel column of radius R = 4 cm; the friction coefficient at 
the walls is fi s — 0.5. 
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FIG. 17: Weight at the bottom of a 2D elastic material, when an overweight equal to the saturation mass is added on top of 
the material, for various Poisson ratio v v : 0.33 (open stars), 0.40 (diamonds), 0.48 (open down triangles), 0.57 (triangles), 0.67 
(open circles) and 0.78 (squares). The friction coefficient at the walls is fi a — 1.0. The data are scaled with the saturation mass 
obtained in a simulation without overweight. 
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FIG. 18: Weight at the bottom of a 3D elastic material, when an overweight equal to the saturation mass is added on top of 
the material, for Poisson ratio v v — 0.35 (triangles), v v — 0.4 (circles) et v v — 0.45 (squares). The friction coefficient at the 
walls is fi s = 0.5. The data are scaled with the saturation mass obtained in a simulation without overweight. 



B. Simulations with an overweight 

In this section, we present simulations of a Janssen experiment, when an overweight corresponding to the saturation 
mass is added on top of the material. 

1. Effect of the Poisson ratio 

On Fig. 1171 and IT81 we plot the apparent mass M a versus the filling mass M/m, rescaled by the saturation mass 
M sa t, respectively in 2D and 3D, for various Poisson ratio. 

We now observe that, contrary to the simulations without any overweight, the results depend strongly on i/ p , i.e. 
on the elastic stress redirection constant K e i. The curves all have the same form: M a increases with Mfm, up to 
a maximum M max , then decreases slowly towards the saturation mass M sat . The relative maximum M max /M sat 
increase with v pi and takes its value for higher Mfm/M sa t. 

The proposed rescaling law, similar to the Janssen one, is not correct. We did not find any rescaling law; we thus 
keep our rescaling as a practical representation of data. 
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FIG. 19: Weight at the bottom of a 2D elastic material, when an overweight equal to the saturation mass is added on top of 
the material, for various friction coefficients /i s : 0.2 (squares), 0.3 (open circles)), 0.4 (triangles), 0.5 (open down triangles), 0.6 
(diamonds), 0.7 (open left triangles), 0.8 (right triangles), 0.9 (open hexagons) and 1.0 (stars). The Poisson ratio is v v = 0.78. 
The data are scaled with the saturation mass obtained in a simulation without overweight. 
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FIG. 20: Weight at the bottom of a 3D elastic material, when an overweight equal to the saturation mass is added on top of 
the material, pour friction coefficient fi s = 0.25 (triangles), fi s = 0.5 (squares) and = 0.8 (circles). The Poisson ratio is 
v v — 0.45. he data are scaled with the saturation mass obtained in a simulation without overweight. 



2. Effect of friction at the walls 

On Fig. ^5] et I2UI we plot the apparent mass M a versus the filling mass Mfm , rescaled by the saturation mass 
M sat , respectively in 2D and 3D, for various friction coefficients at the walls. 

We observe that the relative maximum of the apparent mass M max /M sat increases with the friction coefficient at 
the walls, and takes its value for higher Mfm/M sat . 

This result leads to an apparent paradox. As actually the maximum M max increases when M sa t decreases (i.e. 
when fi s and v v increase), we observe that the more the weight of the grains is screened by the walls, the less the 
weight of the overload is screenedl We will propose an interpretation in the following. 



3. Study of the pressure profile 



We can wonder if all these results remain true for the profile. These features could be due only to the presence of 
a bottom. On Fig. I21I we plot the mean pressure profile when an overweight of mass equal to the saturation mass 
obtained in the simulation of a Janssen experiment is added on top of the material. 
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FIG. 21: Integral of vertical stresses F z vs. depth z, for elastic materials of friction coefficient at the walls fj, s — 0.5, and 
Poisson ratios v v = 0.35 (squares) and v v = 0.45 (circles), when an overweight equal to the saturation mass is added on top of 
the material. We plot F z /F sat vs. K e i[i s z/R where F sat = pgnR s /(2K e ifi s ) is the theoretical saturation value for F z . z/R — 
corresponds to the top of the material. 
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FIG. 22: Integral of vertical stresses F z vs. depth z, for elastic materials of friction coefficient at the walls fj, s — 0.5, and 
Poisson ratios v v = 0.35 (squares) and v v = 0.45 (circles), when an overweight equal to the saturation value F sa t of F z . We 
plot F z /F sa t vs. K e ifj,s z/R where F sa t = pgirR 3 /(2K e i/i s ) is the theoretical saturation value for F z . z/R — corresponds to 
the top of the material. 

We observe the same features as for the simulation of a Janssen experiment. In order to go one step further, we 
now add an overweight which imposes on top of the material a mean pressure equal to the saturation pressure; the 
results are presented on Fig. [23 

Once again, we observe the same features as for the simulation of a Janssen experiment. The Janssen rescaling law 
is still incorrect, and the overshoot effect is more important for higher Poisson ratio and friction coefficient. Moreover, 
on the profile the amplitude of the maximum is more important: for friction coefficient /i s = 0.5 and Poisson ratio 
v v = 0.45, we find that the maximum force on the profile is 1.06 times the saturation force F sat , whereas the maximum 
weighted mass is 1.03 times the saturation mass M sat . 

4- Rescaling with the radius 

We finally verify the rescaling with the radius column for simulations with and without any overweight. This 
rescaling is perfect for both simulations (Fig. I25)> for the extreme radius employed in 16] (R = 1.9 cm et R = 4 cm). 
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FIG. 23: Study of the rescaling with the radius R of the column. We plot A'I a /R 3 vs. Mfm/R 3 for R = 1.9 cm (squares) et 
R = 4 cm (line), for an elastic material of Poisson ratio v v = 0.45 and friction coefficient at the walls /i s = 0.5. 

C. Observable features 

To summarize, several features are experimentally observable if the isotropic elastic theory is valid. 

In a Janssen experiment with free top surface, for a same material, the Janssen constant K, deduced from the 
measured saturation mass M sat , must be higher for higher friction coefficient at the walls. 

In an experiment with M sat overweight on top of the material, the apparent mass M a must increase strongly with 
the hlling mass Mfm, and then decrease slowly towards the saturation mass M sat . The observed maximum M max 
must increase with friction at the walls and with Janssen coefficient K (measured in a Janssen experiment with free 
top surface). 

The saturation mass M sat deduced from the pressure profile must be higher than the one measured at the bottom, 
and a strong decrease in the pressure must be observed near the bottom. This kind of feature would not be observed 
for a hyperbolic theory such as OSL 7]. Moreover, for experiments with an overweight on top of the material, the 
maximum of M a will be higher on the pressure profile. 

D. Comparison with experimental results 

In this section, we compare the elastic theory predictions to the experimental results obtained in [l^. The main 
features have already been presented in 0] 

1. Classical Janssen experiment 

The Janssen experiment data are perfectly fitted by the Janssen model. Therefore, they are also perfectly fitted by 
the elastic theory: we indeed showed in Sec. IIV Al that for weak friction at the wall (/x s = 0.25 as in the experiment) 
the Janssen model and elastic theory predictions cannot be discerned. 

We also showed in Sec. IIV A 31 that when the friction coefficient is increased, the elastic theory predict a sharper 
initial increase of M a with Mfm and an more abrupt saturation. However, this cannot be observed experimentally 
for the small range of friction coefficient used in 0] (from 0.22 to 0.28). In order to test this last prediction, it would 
be necessary to work with higher friction at the walls (/i s ~ 0.5). 

Note that for high packing fractions (y — 0.645 ± 0.005), we obtained Janssen coefficients higher than 1 (K = 
1.2 ± 0.1): this cannot be obtained in the isotropic elastic theory, as K e i = v v j(\ — v p ) < 1. We saw that K e i is less 
than K (Fig. lTJl ll2|) . due to presence of a rigid bottom, but for small friction at the walls, this cannot explain high K 
values: we deduce from Fig. ^]that for /i s = 0.25, K must be less than 1.03; for /x s = 0.8, K could be as high as 1.2. 
We will see in Sec. IIV II how to obtain higher K values in the context of anisotropic elasticity. 
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FIG. 24: a: Apparent mass M a vs. filling mass Mfui, rescaled by the saturation mass M aa t, for loose packing in medium-rough 
columns of 3 diameters (38 mm (squares), 56 mm (circles), 80 mm (triangles)) with an overweight equal to M sa t\ the dotted 
line is the hydrostatic curve, b: Simulation of the experiment of Fig. I24fa ,. for an elastic medium characterized by the same 
saturation mass (Poisson ratio v v — 0.46) and the same friction at the walls (fj, s — 0.25). Inset: maximum mass M max rescaled 
by saturation mass M sa t vs. static coefficient of friction at the walls in experiments made on loose (squares) and dense 



(circles) packing in 38 mm diameter columns, and in simulations for elastic media of Poisson ratios 



0.45 (open squares) 



and v p = 0.49 (open circles); the left vertical scale is used for the experimental data, the right vertical scale is used for the 
simulation data. 



It was observed in [16| that for different preparations, characterized by different packing fractions, the saturation 
mass M sat is lower for higher packing fraction. Although there may be structural differences between the pilings other 
than the packing fraction, this can be interpreted as increase of the Janssen coefficient K with packing fraction. Thus, 
in the context of isotropic elasticity, this leads to an effective increase of the Poisson ratio v v with packing fraction. 
In 0, an effective relation between packing fraction v and Poisson ratio v p was derived: v p ~ 2.i{v — 0.41) with a 
precision of 5%. Note that the largest packing fraction v = 0.645 ± 0.005 would give a Poisson ratio v p = 0.54 ± 0.03 
marginally larger than the limit value of 1/2, as commented above. 



2. Overweight experiment 

The elastic theory gives qualitatively the same behavior as the experimental results (see Fig. I24[) : the apparent 
mass M a increases with the filling mass Mfm , up to a maximum M max , then decreases slowly towards the saturation 
mass M sat . Furthermore, all the features predicted by elasticity can be observed: R 3 rescaling, M max /M sat increase 
with friction at the walls, and M max /M sat increase with packing fraction (experiment) or Poisson ratio (theory). Note 
that M max /M sat increases with packing fraction and Poisson ratio are equivalent because of the effective increase of 
Poisson ratio with packing fraction deduced from the classical Janssen experiment. 

However, there is no quantitative agreement between isotropic elasticity predictions and experiments. Fig. I24f) . 
The elastic curve is very similar to the experimental one, but the experimentally observed maxima are 30 to 40 times 
larger. 
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Rescaled height H/D 

FIG. 25: Mean resistance force to pushing as a function of the height H of the packing scaled by the column diameter D, 
for 1.58 mm steel beads, of packing fraction 62.5%, in a 36 mm duralumin column, for V = 16 fims" 1 (filled squares) and 
V = 100 pirns -1 (open circles). The line and the dotted line are fits by eq. I|28|l . The dashed line is the hydrostatic curve. 



E. An apparent paradox explained by elasticity 

Interestingly, we find in the elastic case the same qualitative phenomenology as in the experiment i.e. the computed 
values of the overshoot M max rescaled by the saturation mass M sat increases both with the friction at the walls and 
the Poisson coefficient (i.e. with the effective Janssen's parameter). This features reads as a paradox: the more the 
weight of the grains is screened by the walls, the less the weight of the overload is screened, but we can now try to 
understand it at least qualitatively. 

If we impose on the top surface of an elastic medium the asymptotic values for displacements and stresses (see eq. 
jll) , (|12fl ). these values then extend to the rest of the column. Thus, with such an overweight, a flat pressure profile 
along depth z a zz {r,z) = — is obtained as in Janssen's theory. Actually, with the overload, the displacement 

imposed experimentally on the surface is almost constant: u z (r) = uo since the overweight is much more rigid than 
the material. Then, as the asymptotic displacement is parabolic, we must have a "transition regime", which is at the 
origin of the overshoot effect. 

There are two limits in which this transition can disappear, i.e. when the Poisson coefficient v or when friction 
at the walls \i s are decreased to zero. Then the parabolic part of the asymptotic displacements (eq. I|12[) ) becomes 
negligible and therefore, the imposed displacement on the surface is close to the asymptotic value. This results in a 
decrease of the overshoot amplitude. 

Basically, we thus recover in the elastic case, the same paradox as the experimental situation. We now understand 
it as a consequence of the boundary condition imposed experimentally by the overweight i.e. an almost constant 
displacement on the surface. 

It would be interesting to put overweights of different Young modulus on top of a granular material in order to see 
if the overshoot amplitude decreases when the overweight Young modulus is decreased. 



F. Pushing experiment 

In this section, we study the elastic theory predictions for the force needed to push an elastic material upwards 
at constant velocity in a column, and compare it to recent experimental observations [LH Il5l | in the case of a slowly 
driven granular material. The main features of this analysis have been presented in (la ] . 

For a vertically pushed granular assembly [T3 . Hll ] at constant velocity, the resistance force F increases very rapidly 
with the packing's height H(see Fig. 

Following the standard Janssen screening picture, this strong resistance to motion is due to the leaning of the 
granular material on the walls (eq. Q) in association with solid friction at the side walls. At the walls, we suppose 
a sliding of the granular material at a velocity Vq (the driving velocity) ; the shearing stress is then 

o-rz(z) = -Vd(Vo)o- rr (z) (27) 

where /Xd(Vb) is the dynamic coefficient of friction between the beads and the cylinder's wall at a velocity Vq. 
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FIG. 26: Comparison of the resistance force to pushing simulated for a homogeneous isotropic elastic medium (squares) of 
Poisson ratio v = 0.45 and Young modulus E = 100 MPa in a duralumin cylinder of radius R = 1.9 cm, with coefficient of 
friction fid = 0.2 at the walls, to the curve obtained with eq. I|28|l with Janssen coefficient K — K e i = — v) = 0.82 (line). 

The force F exerted by the grains on the piston can be derived from equilibrium equations for all slices, thus we 
obtain: 

TT 

F = ggXirR 2 x (exp( — ) - 1) (28) 

A 

where g is the mass density of the granular material, R is the cylinder radius and g the acceleration of gravity. The 
length A = R/2K/id(Vo) is the effective screening length. 

We see on Fig. 1551 that the data are well fitted by eq. ggj. We obtain Kfj, d (V) = 0.140 ± 0.001 at V = 100 (ims 4 
and Knd(V) = 0.146 ±0.001 at V = 100 ^ms" 1 . 

In order to get the isotropic homogeneous elasticity prediction for the pushing experiment, we perform a series of 
numerical simulations using Finite Element Method |24| . The condition a rz = —fi^rr is now imposed everywhere 
at the walls (for the pulling situation, we impose a rz — -\-^a rr ). The cylinder is modelled as a duralumin elastic 
medium. As long as the Young modulus E of the elastic medium is less than 500 MPa, which is usually the case for 
granular media, we find no dependence of the results on E. 

We find no appreciable difference between the elastic prediction (Fig. EEJ) and the curve given by eq. iffify with 
K = K e j. Therefore, regarding the dependence of the stationary state force F on the height of beads, our system 
cannot be distinguished from an elastic medium. 

G. Pushing vs. pulling 

As a check of consistency, we performed the following dynamical experiment in [poT |. First, the granular column is 
pushed upwards in order to mobilize the friction forces downwards and far enough to reach the steady state compacity. 
Starting from this situation, the friction forces at the walls are reversed by moving the piston downwards at a constant 
velocity Vdown = 16 /ims , until a stationary regime is attained. Note that this stationary regime is characterized 
by the same compacity v w 62.5% as in the pushing situation. In Fig. [23 the pushing force F is measured for 
different packing heights H. The fit of experimental results with eq. gives K^^W /ims -1 ) = 0.156 ± 0.002 which 
is 10% larger than ^^^(16 /ims -1 ) extracted from the pushing experiment. This difference, though small, can be 
observed out of uncertainties, and is systematic. It cannot be due to a slight change in compacity v as from relation 
AK/K w 5Ai>/j>, we would expect a 2% variation in compacity between the pushing and the pulling experiment, 
which would be observed; we actually measured Ai) jv = ± 1%. According to Janssen's picture, this would imply 
that vertical stress redirection is more efficient in the downward pulling situation. We believe this is a clear evidence 
of a granular structuring effects but its also shows that this effect is not dominant: it affects only 10% of the average 
mechanical parameter K . 

Note that finite element simulations show that the presence of a rigid bottom implies that the effective Janssen's 
parameter K e ff extracted from Janssen's rescaling for the pulling situation is higher than K e i, whereas for the pushing 
K e ff w K e i (as can be seen on Fig. |23 the fit of the elastic curve with K = K e i is good). Actually, if we adjust 
the elastic predictions for pushing and pulling experiments with an elastic material of Poisson coefficient v p = 0.45, 
eq. (|28|l yields a Janssen's constant K e ff for the pushing which is about 3% lower than K e ff for the pulling. This is 
qualitatively (though not quantitatively) in agreement with the experimental results. 
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Rescaled height H/D 

FIG. 27: Resistance force to pushing (filled squares) and to pulling (open squares) in the steady-sliding regime at V = 16 /ims -1 
as a function of the height H of the packing scaled by the column diameter D, for steel beads in the duralumin column. The 
lines are the fit with eq. 12811 of the resistance force to pushing and its prediction for the pulling situation; the dotted lines are 
the fit with eq. I|28|l of the resistance force to pulling and its prediction for the pushing situation. 



Therefore isotropic elasticity can be a good framework only if we neglect the existence of bulk restructuring effects 
inducing differences in the effective Poisson coefficient of the material between the pulling and the pushing. Note that 
in this case, an isotropic modelling of the granular material is somehow questionable. 



H. Towards anisotropy 

We have seen in the preceding sections that isotropic elasticity reproduces qualitatively all the features observed in 
the experiments performed on granular materials. However, some problems remain: a Janssen constant K of order 
1.2 was observed experimentally whereas elasticity cannot provide Janssen constants higher than 1.03 with the 
same experimental parameters (i.e. /i s < 0.3); moreover, the amplitude of the overshoot predicted by elasticity is 20 
times lower than the one observed experimentally. 

That is why we study here the predictions of the simplest extension of isotropic elasticity: transversely isotropic 
elasticity. 

The stress-strain relations are now: 
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Following the same calculation steps as in 
imposes: 
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Sec. IIII Al . one can show that, in the limit of high depths z, confinement 
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FIG. 28: Apparent mass M a vs. filling mass M/m, rescaled by the saturation mass M sa t, for various transversely isotropic 
elastic media characterized by the same saturation mass M sa t- Squares: G = 0.345 MPa, Ex = E 2 — 1 MPa, v\ = v% = 0.45 
(isotropic medium). Open circles: G = 0.1 MPa, Ex = E 2 = 1 MPa, v x = 1/2 = 0.45. Triangles: G = 1 MPa, Ex = E 2 = 1 
MPa, m. = vi = 0.45. Open down triangles: G = 0.345 MPa, Ex = 2 MPa, E2 = 1 MPa, n = i/ 2 = 0.29. Diamonds: G = 0.345 
MPa, £i = 2 MPa, E 2 = 1 MPa, fi = -0.34, ^2 = 0.55. Open left triangles: G = 0.345 MPa, Ex = 2 MPa, £2 = 1 MPa, 
fi = 0.756, = 0.1. Right triangles: G = 0.345 MPa, Ex = 1 MPa, £ 2 = 2 MPa, i/j = 0.633, v 2 = 0.6. Open stars: G = 0.345 
MPa, Ex = 1 MPa, E 2 = 2 MPa, ^i = 0.815, 1/2 = 0.3. In all these simulations, the friction at the walls is fi s = 0.25. 



i.e. we recover in the anisotropic case a Janssen-like relation between stresses, with a stress redirection constant 

Ei v 2 



The constraints on the elastic parameters are now 



E2I — V1 



v\ > -1 
Ei 

V ^E~ 2 +1 



(36) 

(37) 
(38) 

(39) 



We can see the consequences on K anis through an example. If v\ — 0, then the constraint is v\ < which leads 
to K an i S , < J By adjusting the modulus E\ and E2, one can give any value to K an i S , which is not bounded 

anymore by a maximum value of 1 (the isotropic case). This means that the experimental values found in [Tfij for 
dense packing, i.e. a Janssen constant K ~ 1.2, which could not be reached by the isotropic elastic theory, can be 
understood in the framework of anisotropic elasticity. Note however, that we have now 5 independent parameters 
instead of 2. 

The problem is now that there arc several ways to give K a value by adjusting independently (while satisfying the 
constraints) 4 parameters V\, v 2 , E±, and E%; another independent parameter is G, which may affect the Janssen 
profile shape; thus one has to find a physical justification for choosing these values. 

We tried to vary independently most parameters while keeping K and /i s constant in order to see if for a given 
M sa t, there is a way to obtain the giant overshoot we observe experimentally. 

Several tries are compared with the isotropic case on Fig. EHl As far as we could see, there is only a rather small 
influence of anisotropy on the overshoot amplitude: the deviation from the isotropic value is at the most 15% for 
reasonable values of the modulus. Nevertheless, there are some general tendencies: the higher the shear modulus G, 
the smaller the overshoot amplitude is; the higher V\, the higher the amplitude is, whatever the stiffer direction may 
be. However, the effect is far from sufficient to reproduce experimental results. 



I. A stress induced anisotropy? 



In order to account for the height of the experimentally observed overshoot, we propose a toy model based on the 
idea of stress-induced anisotropy. It is possible that the overweight induces locally a change in the structure. We 
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would expect an increase of the number of contacts i.e. of the young modulus E2 in the vertical direction, if the 
granular material can be modelled as an effective elastic medium. From the relation K anis = ^ ^ 2 , we then expect 
the Janssen coefficient to be lower near the overweight than in the bulk. 

Let the Janssen constant be K\ near the overweight, down to a depth Hi, and K% > K\ for any depth z > Hi. The 
saturation mass, which is the overweight mass, was measured with an homogeneous material of Janssen constant K2 
and is thus M sat2 = , Then, it underestimates the saturation mass M satl — of the material layer near 

the overweight: the measured mass M a then first increases with the filling mass Mfm and would naturally tend to 
a higher value M sa t 1 . But for a depth Hi the material's structure is no more affected by the overweight and is now 
characterized by a Janssen constant K^: the mass imposed on the material 2 by the material 1 is then higher than 
its saturation mass M sa t 2 , and the measured mass M a has to decrease with the filling mass Mfm in order to reach 
its saturation value M sa t 2 ■ 

This simple idea can be easily formalized for an ideal Janssen material. The equilibrium equation on the vertical 
stress a zz (z) in material i characterized by Janssen constant Ki reads 

da zz 2Kip s 

+ — 5— °zz = -P9 (40) 
oz R 

In the material 1, for < z < Hi, we solve this equation for the mass M(z) = <j zz (z)/(iTR 2 g) weighted at depth z, 
with the boundary condition M(z = Q) = M sa t 2 and get 

M(z) = M satl (l-exp(-2K 1 v s ^yj + M sat2 exp(-2ifi/z s -|) (41) 

In the material 2, for Hi < z we solve this equation with the boundary condition M(z = Hi) = M satl (l — 
exp(-2Kifi s ^)) +M sat2 exp(-2ir 1 /i 5 ^) and get 



(42) 



M(z) = M sat2 + (M satl - M sat2 )(l - exp(-2^ lMs i|)^ exp(-2if 2M , (z 
These equations can be rewritten using M a /M sat2 and Mfni/M sat2 as variables: 

/ rr M a M sa tl {, ( M fiU M sata \\ ( Mfm M sat2 \ . . 

< z < H 1 — > — = — 1 - cxp -—r — — + exp — — f- — — 43 
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„ . M a f M satl ,f ( Hi M sat2 ^\ / . M fiU Hi \ 

The two independent variables one can adjust are jf 2 ^ 1 (which is the Janssen constant ratio K2/K1), and the ratio 

H1/X2 = Mfui(Hi)/M sa t 2 where A2 = pirR/ '(2K2H S ) is the Janssen screening length in the medium 2. 

On a M a /M sat2 = f (M f m / M sat2 ) plot, the parameter H1/X2 is the X-axis value Mfui/M sat2 at which the weighted 

mass starts to decrease. The slope at the origin is 1 — M s "| 2 . Typically (see Fig. l2*$land Fig. l37)|) . the experimental 

slopes are of order 0.5 (it would be 1 for an hydrostatic pressure), i.e. the Janssen coefficient Ki has to be about 
2 times K2, in the anisotropic elasticity framework, it would mean that the young modulus in the vertical direction 
is doubled by the overweight. The order of the extension H1/X2 of the induced anisotropy is of order 0.4. The 
experimental data for dense and loose packing are compared on Fig. 1291 and Fig. 1201 with the model. It appears to 
reproduce correctly the data; however, the parameters cannot be determined accurately as they are obtained from the 
very beginning of the data (the slope at the origin and the X-axis value of the maximum) which cannot be measured 
with a high precision. Therefore, the parameters presented in these figures are just given as examples. 



V. CONCLUSION 



In conclusion, we performed an extensive study of the Janssen's column problem. It is a mixture of numerical 
studies both in 2D and in 3D, in the case of frictional boundaries. The aim was to test thoroughly the classical and 
celebrated Janssen's analysis in the context of an effective homogeneous elastic material and provide some meaning 
to the effective Janssen's constant of stress redirection at it can be obtained experimentally. Note that here we do not 
make any assumption on a plastic threshold in the bulk as it is usually considered to provide bounds on the Janssen's 
constant (active and passive limits). Also this analysis was performed in the context of an extensive experimental 
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FIG. 29: Apparent mass M a vs. filling mass Mfui, rescaled by the saturation mass Msat, for dense packing in medium-rough 
columns of 3 diameters (38 mm (squares), 56 mm (circles), 80 mm (triangles)) with an overweight equal to Msat- The data are 
compared with our inhomogeneous Janssen model with K1/K2 = 2.4, and H1/X2 = 0.45. 
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FIG. 30: Apparent mass M a vs. filling mass Mfm, rescaled by the saturation mass M 3a t, for loose packing in medium-rough 
columns of 3 diameters (38 mm (squares), 56 mm (circles), 80 mm (triangles)) with an overweight equal to Msat- The data are 
compared with our inhomogeneous Janssen model with K1/K2 = 2.8, and H\/\2 = 0.3. 



work where preparation was varied and were special care about friction at the wall was taken to be able to establish 
a precise comparison with theoretical modelling. Interestingly, we find that the Janssen's approach is fully valid in 
the limit of low friction coefficients between the grains and the wall (up to a moderate value of 0.5). It means that an 
exponential saturation curve of the average normal stress at the bottom is an excellent approximation, which defines 
precisely an effective Janssen's constant solely dependent on the Poisson ratio. We also derive ion the limit of an 
infinite column a value for an elastic Janssen's constant (K e i = v in 2D and K e i = v/(l — v) in 3D). For a finite size 
column, the presence of a bottom diminish the average vertical stress such as to yield an effective Janssen's constant 
if the bottom normal stress is measured with a value directly related to the elastic constant K e i. Therefore, in this 
context experiment data can be matched with isotropic elasticity if the packing fraction representing the preparation 
can be associated with an effective Poisson ratio. 

Consequently, there is a need to provide a more strained test to the isotropic elastic theory. This was done 
experimentally by imposing on the top of the column an overweight equal to the saturation stress |l6j and here we 
propose the same test in the same conditions for an elastic material. The numerical simulations show that stresses 
at the bottom also exhibits an overshoot when the column depth is increased. The relative value of the overshoot is 
at most 7% . We propose a scaling relation of the overshoot amplitude when the wall fiction is varied. This result 
contrasts with a Janssen's analysis where a fiat profile should be observed. This overshoot effect is related to the 
deformation of the elastic medium below the overweight which sets a length where the deformation profile is not 
parabolic any more as in the rest (more like a flat profile as it is imposed by the overweight boundary condition). 
But at the quantitative level, the overshoot effect found experimentally has an amplitude about 20 to 30 times larger! 
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This lack of quantitative agreement opens new questions on the modification on the influence of the medium due 
to the overload. This is the reason why we push further the investigation in the context of anisotropic elasticity. 
We consider an orthotropic elastic medium with the main stiff direction along the vertical. In this case the effective 
Janssen's redirection coefficient can be changed according to the stiffness ratio but the overshoot test does not produce 
an overshoot value significantly larger than the isotropic situation. 

Finally, we propose a qualitative model based on an extension of the Janssen's approach where we assume that the 
overload has changed the medium within a given depth such as to yield a smaller Janssen's constant. This would 
be consistent with the onset of strain induced anisotropy producing a stiffer medium in the vertical direction. The 
agreement of this simple model (with two fitting parameters) is satisfactory but more importantly it raises interesting 
questions and calls for new experimental work. In this frame of mind it would be very interesting to see differences 
of the stress saturation curve for two media prepared i) with a regular rain like pouring as before and ii) a rain 
like deposition process but where an overweight is imposed above each deposition step (a deposition step being of a 
height much smaller than the final height). In this case the stress induced texture changes could provide a rational 
explanation to the extension of the Janssen's model that fits correctly the overweight experiments. 

We thank Profs. R.P. Bchringer and J. Socolar for many fruitful discussions. 
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